Canonical-grandcanonical ensemble in-equivalence in Fermi systems? 
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I discuss the effects of fermionic condensation in systems of constant density of states. I show that 
the condensation leads to a correction of the chemical potential and of the Fermi distribution in 
canonical Fermi systems at low temperatures. This implies that the canonical and grandcanonical 
ensembles are not equivalent even for Fermi systems. 



INTRODUCTION 

Let a system S be in contact with a heat and particle 
reservoir, TZ. The microstates of the system will be de- 
noted by TOi. If I assume that the system is ergodic and 
all the microstates corresponding to any fixed E and N 
are equally probable, then the probability associated to 
any state, say p{mE,N) satisfies 



p{mE,N) oc e^(^-''^) 



(1) 



where (3 = l/^ksT), T is the temperature, and /i is the 
chemical potential of the reservoir. 

The same probability distribution may be obtained if 
we take as starting point information theory [3| . The in- 
formation we start with is that the average energy and 
particle number are E and A'^, respectively. Then, the 
probability distribution (Q) is the least biased estimate 
possible on the given information If the system is 
ergodic or not, is of no importance from this point of 
view. For a macroscopic system the maximum of the 
probability distribution is sharply peaked around the av- 
erage values, {N) and (E), so both, particle number and 
internal energy, have well defined values. 

Now let us assume that we can calculate (and measure) 
a macroscopic parameter, which I shall denote by X. The 
probability distribution over the microstates is given by 
Eq. ^ and from this one can calculate the probability 
distribution over the parameter X: 



X{mi)=X 

Px{X) - Pi"^^ 



(2) 



If X is well defined, then Px{X) should have a sharp peak 
aX X = (X). If Px{X) has two maxims, then the system 
undergoes a phase transition. In each of the phases, one 
maximum dominates and this fixes the value of X for 
that phase 0,0]. 

For a system of fermions a parameter which is sur- 
prisingly interesting to analyze is the number of particles 
that occupy completely an energy interval (no holes left 
in this interval) jStarting at the bottom of the single par- 
ticle spectrum H, ^ |^ - let me call this parameter Nq 
and the energy interval [0, eo]. In Q I gave an example 
of an interacting system for which the probability distri- 
bution P/Vo(-^o) forms, below a certain temperature, two 



competing maxims. One maximum, which exists for any 
T > 0, is located at A'o = and the other appears at 
finite Nq. At transition temperature the maximum cen- 
tered at Ao > equals the maximum existent at Aq = 
and a first order phase transition occurs. Above transi- 
tion temperature, since P;vo(A^o) is maximum at A'o — 0; 
(Ao) (^ 0) is microscopical. Below transition temper- 
ature P{No) is maximum at Nq > 0, so (A^o) > is a 
macroscopic quantity. Due to the interaction, an energy 
gap is formed between the degenerate A'o particles and 
the rest of the particles. 

The same parameter may be analyzed for a system of 
ideal fermions. Assume that the density of states (DOS) 
has the general form cr(e) = Ce*, where C and s are 
constants. Now I require that Aq particles form a de- 
generate subsystem on the first A^o energy levels, and the 
first hole in the spectrum appears at energy eg (or energy 
level A'o -I- 1 - see Ref. II for details). Using again Eq. 
(O, I calculate PNg{No) or P^gieo). If Z is the partition 
function of the system, and Zf^^ is the number of con- 
figurations with the first hole appearing at A'o + 1, then 
Pno{No) = ZnJZ and P,„(eo) = fT(eo) • [Ao(eo)]. 
Since Z is a, constant, the extrema of P are found by 
solving dZj\jg/dNQ = 0, or dZ^^/deo — 0. Even more 
convenient is to work with logZ^rQ, which is 



log-ZjVo 



s - 



s+2 



+C 



2 

de e" log 



1 



(3) 



Since dZj\jg/dNo = {dZj^g/ de^) ■ {deo/dNQ) = 
{dZj\jg/deo) ■ cr^^(eo) = implies dZxo/deo = 0, I cal- 
culate H 



dlogZ 



No 



den 



{log 










C4l 



= 0. 



(4) 

If s > 0, logP has one and only one maximum at 
eo,A^Q > 0, so for any macroscopic systems (i.e. large 
enough C) there will be a degenerate subsystem on the 
lowest energy levels at any temperature. If s = (e.g. 
particles in a two dimensional flat potential or in a one- 
dimensional harmonic potential) there is a transition 
temperature, Tc,2D, below which the maximum of P{Nq) 
moves from Aq = to Aq > 0, i.e. a degenerate gas 



2 



forms. The degenerate gas may be put in correspon- 
dence with the the Bose-Einstein condensate in a gas of 
bosons with similar spectrum H, Q |E 1^ and for simplic- 
ity I shall call it the Fermi condensate or the degenerate 
subsystem. The most interesting case seems to be s < 0, 
when logP has either only one maximum, at A'o = 0, or 
two maxima, at iVo = and A^o > 0. As the temperature 
decreases, the second maximum increases and becomes 
bigger than the maximum at A^o = 0- 

The rest of the article is organized as follows: in Sec- 
tion 2 I discuss in some detail the probability distribution 
P/Vq(A^o) for the two-dimensional Fermi gas and calcu- 
late the relevant parameters of the Fermi condensate. In 
Section 3 I use in the standard way the grandcanonical 
formalism to calculate the thermodynamic quantities. In 
section 4 I compare the entropy obtained in section 3 to 
the entropy of the thermodynamically equivalent Bose 
gas and point out the miss-fit from the low temperatures 
range. The solution is suggested in section 5: in the 
canonical or microcanonical ensemble the condensate re- 
gion should be considered separately and the Fermi dis- 
tribution applies only in the thermally active layer, which 
is the energy interval above the condensate. The last sec- 
tion is reserved for conclusions. 



FERMI CONDENSATION AT CONSTANT 
DENSITY OF STATES 

For constant a, Eq. Q has a solution if and only if 
log[l -I- e-i^''] < (aksT)-'^. Therefore, as in jj], I define 
the condensation temperature Tc^f by the equation 



log[l 



-/3c,fM1 — 



= [crkBTcFY 



(5) 



If I assume that in the temperature range of interest 
aksT :$> 1 (i. e. I consider macroscopic systems), then 
also Pfj, ^ 1, and the equation above may be written as 



(6) 



In Eq. © I used the approximation N w a^, valid at 
low temperatures (see further, Eq. O. For T < Tc^f, 
1 > log[l + e'^^'^o"^)] « e'^^'^o^^) and the solution of Eq. 
Q may be approximated by 



eo,r,iax = - /csT log [crfcsT] 



(7) 



The particle number that is associated to eo,ma£c is 
No^max = ea.rnaxO- = a ji-aksT loglaksT]. The distribu- 
tion is not symmetric and (iVo) < No^max- In the low 
temperature limit, (A'o) converges to Ng^max, but gen- 
eral, closed analytical expressions for (A'o) seem difficult 
to find. 

Let me denote (eo) = {NQ)/a. In the next sections 
we shall see that in a canonical finite system we have 
to consider that all the energy levels from to (eo) ^ 



^o,max are occupied, and the Fermi distribution applies 
only to the energy levels from (eg) up-wards. For this 
reason, as mentioned in the Introduction, I shall say that 
the particles in the energy interval above (eq) form the 
thermally active layer. Obviously, both (cq) and (Nq) 
are subject to fluctuations. By comparing Eq. {T)) or 
Q with Eq. lO, we observe that T is the condensation 
temperature for the gas in the thermally active layer. The 
total number of particles in the system is then calculated 
as 



N 



{No) 



a de 



{eo)+<T-i & 



1 



= cr/i — 1 



+akBT log 



l4.e^((^o)+<T-i-A0 
akBTe^^^"'^+''''-^''^ Ka^i + e^ 

1 W cr^ + /3((eo) - eo,niax + d 

/3((eo) -eo,ma.) + e'^^"' 



CT/i — 1 

/3((eo)-co.™ax-|-o-i) 



UEF ; 



(8) 



where e^? is the Fermi energy and I used the fact that 
crfesTe'^^^o-""--'^) = 1 (Eq. |i|or[7|). 

In the grandcanonical ensemble, if we do not take 
into account the degenerate subsystem, the total particle 
number is 



iV = cr/i + crfcsT log [1 + e"^^] 



(9) 



The Fermi energy, denoted in this case by , is given by 
the equation = cr^F, so for ^ 1 we have 



eF ~ l-J- 



keTe-^'' 



(10) 



GRANDCANONICAL FERMI GAS 



Under grandcanonical conditions, the entropy and in- 
ternal energy of a Fermi gas have the expressions: 

S = -aklT [2L^2i-e^^') + f3fi\og{l + e^^')] (11) 
and U = -{kBTfaLi2i-e^'') , (12) 

respectively. If we assume that different ensembles are 
equivalent, then we can use Eqs. © and (|12|l to express 
fi and T in terms of N and U and plug these expres- 
sions into S. The result should be the same. Vice-versa, 
calculating T = (dS/dU)-^ and -^/T = {dS/dN)-\ 
one should obtain the grandcanonical temperature and 
chemical potential. 

To illustrate this change of variables, I shall express 
S (Eq. Ill|l in terms of U and N in the limit of low 
temperatures. In this limit ^ 1 and I can neglect 
Q-0t^ from Eq. (|5J|, retaining 



f3fi. 

On the other hand, using the expansion 



3(/?/l)2 



(13) 



(14) 
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FIG. 1: Surface L the entropy of a Fermi system, s = S/kso, 
as a function of A'^ = N/a and u = U/a. The range on the 
vertical axis is from to 7r-y/2/3 (see Eq. I16II . Surface II is 
obtained from surface I by substracting from the energy (7, 
the ground state energy of the system, Ug.s.{N) = /2a. 



and Eq. (|13|l . U may be approximated as 

2' 



U 



2^ 



N 



(15) 



Eliminating fc^T and (3fi from Eqs. H13() and (|15|l plug 
ging it into Hll|l I get 



S 2 fU 



(16) 



The result (|16|l is plotted as surface I in Fig. ^ 

Since Uo{N) = N'^/2a is the zero temperature en- 
ergy of the Fermi gas, let me denote the excitation en- 
ergy hy Ub = U - Uo{N) and define S'b(C/b,7V) = 
S{Ub + Uo{N),N). A gas that has the internal energy 
Ub and entropy Sb{Ub,N) as defined above, is called 
thermodynamically equivalent with the original Fermi gas 
(equivalence classes were defined in Q). The chemical 
potential of the gas B is related to the chemical poten- 
tial of the Fermi gas hy ^b = — dUo/dN = ^ — ep- In 
the low temperature limit 



Sb 2 Ub 

kB<J V 3 (T 



(17) 



which is plotted as surface II in Fig. ^ Another expres- 
sion for U B may be obtained by using Landen's relation, 
U2{-y)+Li2[v/{v + l)] = ~\\og\l + y) 

UB = {kBTYLi2[{l + e-P''r^]. (18) 

Applying Landen's relation to Eq. (|11|) . I get 

S = _afc|T[2Li2(-e^^)+/3Mlog(l + e'3^)] .(19) 



THE EQUIVALENT BOSE GAS 

It is well known that all the ideal gases of the same 
constant DOS and anyexclusion statistics fall in the same 
equivalence class 1^ IH S S ^i 0| ■ Since the ground 
state energy of a Bose gas (energy at zero temperature) is 
zero, Ub defined in the previous section is nothing but the 
energy of the Bose gas and the surface II in Fig. ^is the 
entropy of the Bose gas. From the bosonic perspective. 



Ub^ 



ade 



p/3(e-AiB) _ I 



a{kBTfLi2{e^^n. (20) 



and from Eqs. ^ and UHl I obtain e"'^^^ = I + e" 

In principle a macroscopic Bose gas of constant DOS 
does not condense, but in a finite system the ground state 
becomes macroscopically populated below a "condensa- 
tion" temperature, r^s- The transition is not sharp, so 
T^.B may only be estimated. Grossmann and Holthaus 
used the equation 



N = akBT,^BlogN 



(21) 



as definition for r^s- Comparing T^.b with Tc.f, one 
can easily see that Tc^f > Tc,b, but they differ only by a 
factor Tf/Tb = l+\og[\ogiakBTc^F)]/ \og{akBTc,F) > 1. 

Let me now examine closer the equivalence between 
the Bose and Fermi gasses below Tc^f- At these tem- 
peratures, the ground state population of the Bose gas 



1 



1 



akBT,,F, (22) 



which is of the order of N. In such a case the population 
of the ground state should be considered separately in 
all the equations. The total particle number in the Bose 
system should then be written as 



N 



No 



a de 



-1 e 



/3(e-MB) _ 1 



No 



-CrfcsTlog 1 



(23) 



Much below Tc^f, cr/is ^ 1 and Eq. 
proximated by 



may be ap- 



N-No~akBT\og{akBT), 



(24) 



which is equivalent to Eq. ((T)). 

Now a natural question arises: does the thermody- 
namic equivalence still holds at temperatures below con- 
densation? For example at a temperature T = rc,s/2, 
from Eq. ^ I get « {2an/N)logN w 21ogiV, 
which, if plugged into the expression H22|) for Nq, gives 



No « e2i°s^ ^N^:$>Nl 



(25) 



One way to interpret this is that the thermodynamic 
equivalence breaks down at temperatures around and be- 
low Tc.F, i-e. when the population of the ground state 
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becomes macroscopic. Therefore, below T^^f the trans- 
formation U Ub and S — > Sb, presented in the pre- 
vious section, does not lead to the bosonic entropy. But 
this cannot be true, since it was shown that there exist 
a one-to-one correspondence between the microscopical 
configurations of a Bose and a Fermi system, of identical 
constant DOS and the same excitation energy j3, 
Therefore, the canonical entropies and partition func- 
tions ought to be the same. In this case, what do we 
have to tailor? The answer is: the canonical partition 
function of the Fermi gas. 



CANONICAL FERMI GAS 

By the exclusion statistics transformation (EST) 
method 0, 0] every distribution of fermions along the 
single particle energy axis is transformed into a distri- 
bution of bosons. The A^'o fermions that form the Fermi 
condensate are mapped onto bosons on the ground state 
(the Bose condensate). For a system of constant DOS 
the rest of the particles have a Fermi distribution and 
are mapped onto a Bose distribution of particles. This 
transformation have been made explicitly in Section V.A 
of Ref. The bosonic chemical potential is given by 
Eq. (113, with No = (No) = {e-l^f^ - l)'^. This fixes 
also the fermionic chemical potential to /i = ep + ^Jt-B 
and the number of fermions in the condensate to (A^o)- 
As a consequence of this, in the limit of low tempera- 
tures /3(eF — /^) does not have the expected asymptotic 
behavior (i{ep — /i) « e"'^'^^, but instead, from Eq. 
I get 



PiiB ^[N- akBT\og{akBT)]-^ 
1 crkBT\og{(jkBT) 
N^ iV2 ■ 



(26) 



In addition, Jaynes' theory - according to which the 
probability associated to each microscopic configuration 
should have the form ^ - seems to contain intrinsic 
contradictions. The distribution iQJ represents the least 
biased estimate given that the average number of parti- 
cles in the system is N and the average energy is U. On 
the other hand, the bosonic distribution of the equivalent 
Bose gas is the least biased estimate, given the average 
excitation energy Ub and particle number N. As shown 
above, the two distributions do not map onto each-other 
over the whole spectrum. In Ref. |5| it is proven that 
the situation is even more dramatic when the DOS is not 
constant. 



CONCLUSIONS 

In conclusion I discussed the effect of fermionic con- 
densation - which is the apparition of a degenerate sub- 



system at the bottom of the single particle spectrum - 
in a system of constant density of states. This leads to 
a correction in the calculation of the chemical potential 
in the canonical ensemble at low temperatures (see Eq. 
I26|l . which further implies that the canonical and grand- 
canonical ensembles are not equivalent at low tempera- 
tures even for Fermi systems. 



By applying the exclusion statistics transformation to 
the Fermi system, one obtains a (thermodynamically 
equivalent (3|) Bose system. If the Fermi system is 
condensed, the degenerate subsystem is mapped onto 
the Bose-Einstein condensate - from where the name of 
Fermi condensate is derived. The condensate fluctuates 
and the value of these fluctuations are simply given by 
the condensate fluctuations of the equivalent Bose sys- 
tem. For the calculation of Bose ground state fluctua- 
tions in canonical end microcanonical ensembles I refer 
the reader to the articles of Holthaus et al. (e. g. P , 
and citations therein) and Tran et al. (e.g. and 
citations therein). 



Non-equivalence between grandcanonical Bose and 
Fermi gases have also been observed very recently by Pat- 
ton et al. |l6l| in computer simulations of small systems. 



[1 
[2: 

[s: 
[4: 
[s: 
[6: 
[7: 

[8 
[9 
[10: 

[11 
[12 

[13 
[14 

[15 
[16 



E. T. Jaynes, Phys. Rev. 106, 620 (1957); Phys. Rev. 
108, 171 (1957). 

K. Lee, Phys. Rev. E 53, 6558 (1996); J. Lee and K.-C. 
Lee, Phys. Rev. E 62, 4558 (2000). 

D. V. Anghel, J. Phys. A: Math. Gen. 35, 7255-7267 
(2002), cond-mat/0105089 

D. V. Anghel, J. Phys. A: Math. Gen. 36, L577(2003); 
cond-mat/0310248 

D. V. Anghel, submitted to J. Math. Phys.; 
cond-mat/0310377 

S. Grossmann and M. Holthaus, Phys. Rev. E 54, 3495 
(1996). 

B. DeMarco, S. B. Papp, and D. S. Jin, Phys. Rev. Lett. 
86, 5409 (2001). 

R. M. May, Phys. Rev. 135, A1515 (1964). 

M. H. Lee, Phys. Rev. E 55, 1518-1520 (1997). 

M. Apostol, Phys. Rev. E 56, 4854 (1997). 

M. Crescimanno, and A. S. Landsberg, Phys. Rev. A 63, 

35601 (2001); cond-mat/0003020 

L. Lewin, Dilogarithms and associated functions (McDon- 
ald, London, 1958). 

M. H. Lee, J. Math. Phys. 36, 1217 (1995). 

M. Hohhaus and E. Kalinowski, Ann. Phys. (N.Y.) 276, 

321 (1999). 

M. N. Tran, J. Phys. A, Math. Gen. 36, 961 (2003). 
K. R. Patton, M . R. GeUer, and M. P. Blencowe, 
|cond-mat/0406285| 



